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Abstract 

Numerical evidence of directed transport driven by symmetric Levy noise in time-independent 
ratchet potentials in the absence of an external tilting force is presented. The results are based 
on the numerical solution of the fractional Fokker-Planck equation in a periodic potential and 
the corresponding Langevin equation with Levy noise. The Levy noise drives the system out of 
thermodynamic equilibrium and an up-hill net current is generated. For small values of the noise 
intensity there is an optimal value of the Levy noise index yielding the maximum current. The 
direction and magnitude of the current can be manipulated by changing the Levy noise asymmetry 
and the potential asymmetry. 
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Fluctuation-driven transport is a topic of considerable theoretical and practical impor- 
tance. Its simplest realization is the Brownian motion in which fluctuations give rise to 
diffusion. In the presence of external forces, fluctuations can have a highly nontrivial effect. 
When inertia is neglected, this problem can be described with an stochastic, overdamped 
Langevin equations of the form 

x = -d,v+F+m. (1) 

where V is the potential, F is an external tilting force, and ^ represents the fluctualtions. 
Of particular interest is the case of periodic potentials, V{x + L) = V{x), lacking reflection 
symmetry. These potentials, for which there is not an xq such that V{—x) = V{x + xq), are 
known as ratchets (see Fig. 1). The spatial symmetry breaking implies a preferred direction 
along which the restoring force on a particle at a potential minima is weaker. Because of this, 
it is tempting to think that in the presence of fluctuations, particles will tend to drift towards 
the weak part of the potential giving rise to a "down-hill" unidirectional net transport. How- 
ever, it can be shown that, consistent with the second law of thermodynamics, such directed 
transport can not take place in the case of fluctuations in thermodynamic equilibrium (i.e., 
for ^ corresponding to thermal Gaussian noise), F = 0, and time independent potentials, 
V = V{x). An appealing discussion that illustrates this issue in a ratchet and pawl me- 
chanical device can be found in Ref.[l|. Nevertheless, in the presence of non-equilibrium 
perturbations, fluctuations can give rise to net directed transport in ratchet potentials. The 
realization and exploitation of this idea has generated a great amount of interest in the non- 
equilibrium statistical mechanics of "Brownian motors" spanning more than a decade, see 
for example Ref. and references therein. A trivial way to generate directed transport is 
by imposing biased perturbations, e.g. a constant F. However, what is highly nontrivial is 
the appearance of directed currents in the presence of symmetric, unbiased, i.e. zero-mean, 
non-equilibrium perturbations. Referring to Eq. ([1]), non-equilibrium perturbations can be 
introduced in F, V and ^, and ratchet models can be classifled depending on how this is 
done. For example, in "pulsating ratchets" V = V{x,f{t)), in "tilting ratchets" a time 
dependent zero-mean tilting force F = F{t) is assumed, and in "temperature ratchets" the 
noise strength T in < ^(t)^(s) >= 2r]kBT{t)6{t — s) is a function of time. A less studied 
class of ratchets, which is the main focus of the present work, incorporate non-equilibrium 
fluctuations by assuming non-Gaussian noises. Some examples include ratchet transport 
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driven by symmetric Poissonian white shot noise with exponential distributed amphtudes 
[sl, and the use of colored non- Gaussian noise 1^. 

In this letter we consider directed transport in ratchet potentials driven by Levy noise 
which gives rise to anomalously large particle displacements known as Levy flights. A 
Levy flight process is a non-Gaussian, non-stationary random process whose increments are 
independent and distributed according to an a-stable Levy probability distribution function 
B. a-.UUe distributions plav a very .mportant role because, according to the generalized 
central limit theorem [g] , they are the attractors of distributions of sums of random variables 
with divergent second moments . Also, these distributions exhibit slowly decaying tails 
describing random processes with large events, and their self-similarity properties make 
them useful in the description of scale free transport and random fractal processes. These 
properties explain the ubiquitousness of Levy statistic in many areas of science, engineering, 
and economics [3]. The role of Levy flights in confining potentials has been addressed in 
the literature before. Examples include the study of the decay properties of the probability 
distribution function, the study of unimodal- multimodal bifurcations during relaxation [Sj, 
and the barrier crossing Kramers problem j^. However, despite the widespread recognition 
of the importance of Levy processes, the effect of Levy noise in ratchet potentials has not 
been addressed. 

Beyond its intrinsic theoretical interest, one of our motivations for the study presented 
here is the problem of non-diffusive transport in plasmas. In a recent letter [l^ a ratchet-type 
transport model was proposed to describe non-diffusive impurity transport in magnetically 
confined fusion plasmas. In that work it was assumed that the fluctuating electrostatic 
potential is a stationary, homogeneous Gaussian process, and thus in order to introduce a 
non-equilibrium fluctuation an additional time variation of the two-point Eulerian correla- 
tion function had to be assumed. Although these assumptions are reasonable and physical 
realizable, recent experiments and numerical simulations have shown that there are cases in 
which fluctuations in turbulent plasmas exhibit Levy statistics. For example. Levy statistics 



larators and 
Accordingly, 



has been observed in electrostatic edge turbulence in tokamaks and ste 
in numerical simulations of pressure-gradient driven plasma turbulence 
an open problem of interest is to explore the role of Levy noise in ratchet-type transport 
of impurities. Our results provide indirect evidence that in the presence of non-Gaussian 
Levy fluctuations, a pinch effect might be present in the case of static potentials. The ideas 
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presented here might also have an impact on ratchet transport in biological systems and 
condense matter in general where there is evidence of non-Gaussian fluctuation phenomena. 

Our methodology consists of two complimentary approaches based on the Langevin equa- 
tion ([T]) driven by Levy noise, and in the solution of the corresponding space-Fractional 
Fokker Planck (FFP) equation. For a discusion on the equivalence of the two formulations 
see Ref. {l3|]. Although the numerical methods for the solution of fractional differential 
equations is a rapidly developing fleld, the integration of the FFP equation with fractional 
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derivatives in space and periodic potentials has not been treated before. In Ref. 
numerical study was presented of the FFP equation in a tilted periodic potential in the 
presence of subdiffusion. The key difference between that work and our contribution is that 
while Ref. [l^ deals with subdiffusive processes caused by non-Markovian memory effects 
modeled with fractional derivatives in time, here we consider superdiffusive processes caused 
by Levy flights modeled with fractional derivatives in space. The ratchet potential is given 
by 

{1 — cos [vrx/ail , < x < Oi 

L / iJ , - 1 

1 + cos [7r(x — ai)/a-2\ , ai < x < L , 

where Vq is the amplitude, L = ai + a2 is the period, V{x + L) = V{x), and A = (01—02)/!/ 
is the asymmetry parameter. In all the calculations presented here, Vo = L = 1. Compared 
with the potential V = Vo[sin(27rx/L) + 0.25 sin(47rx/L)] typically used in the literature, 
e.g. Ref j^, Eq. offers the advantage of an easier control of the spatial asymmetry. As 
shown in Fig. 1, for A = —0.274 both potentials coincide. In the Langevin approach we 
assume that the noise satisfles L{At) = ^{t')dt' where L is a Levy flight process with 
characteristic function, pL{k, At) = e*^^pi(x, At)dx given by [s] 



pLik, At) = exp {-Atx|A;r [1 - isign(A;)/3tan(7ra/2)]} , (3) 

where a is the Levy index, j3 is the skewness parameter, and x is the Levy noise intensity. 
Unless explicitly indicated, we will consider x = 1/2. Since for the ratchet problem the mean 
(flrst moment) of the displacement must be flnite, the range for admissible a is 1 < a < 2, 
whereas —1 < 3 < 1. The time step used in Langevin simulation was At = 10"'^. Details of 
the numerical scheme can be found in Ref. [15[. For a = 2, Eq. ([3]) reduces to a Gaussian 
corresponding to a diffusive (Brownian) process. As expected, in this case the problem 
reduces to the Smoluchowski-Feynman ratchet which, as the insert in Fig. 1 shows, only 
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exhibits a current in the trivial case of a biased constant force F. 

One of the main obstacles in the numerical integration of equations containing frac- 
tional derivatives in space is the fact that the left, a-D", and the right, x-D^, Riemann- 
Liouville fractional derivaties are in general singular at the x = a and the x = b boundaries 
respectively p^]. Here we circumvent this problem by regularazing the fractional space 
derivatives using the Caputo prescription, ^D^P =a [P{x) — P{a) — P'{a){x — a)] and 
ID^P =x [P{x) - P{b) + P'{b){b - x)], and write the FFP equation as 

dtP = [Pd^V^jf] + X [1:d: + r^D^] P , (4) 

where K// = V{x) — Fx, and 

"^^•^"r(2-«) Ja {x-yT-'^^' 

for 1 < a < 2, with P" = dyP{y,t). The expression for %D'^P follows from Eq. ([5]) by 
changing the integration limits to and interchanging x and y in the denominator. The 
factors / = —(1 + /3)/[2 cos(a;7r/2)] and r = —(1 — /3)/[2 cos(Q;7r/2)] determine the relative 
weight of the left and the right fractional derivatives. The order of the fractional operators 
a, the asymmetry parameter /5, and the diffusivity x ^'^^ given by the index, the skewness, 
and the scale factor of the characteristic function of the Langevin Levy noise in Eq. (I3l). 
For the numerical integration of Eq. (jl]) we used a finite difference scheme based on the 
Grunwald-Letnikov representation of the regularized fractional operators. Further details 
of the method can be found in Ref. 17|]. The integration domain, x G (a, &) = (—10, 10), 
covered twenty periods of the potential. We used a grid size Ax = 0.02 and an integration 
time step At = 1.8 x 10"'^. The boundary conditions were P{a) = P{b) = 0, and the initial 
condition was a higly peaked normalized Gaussian distribution centered at x = 0. 

Figure 2 shows typical trajectories obtained from the Langevin equation for different Levy 
noise indices a and degrees of potential asymmetry, A. Panel (a) illustrates the absence of 
directed current in the Gaussian case, whereas panel (b) shows a positive bias in the presence 
of symmetric Levy flights. As expected, the Levy noise gives rise to large particle jumps 
encompassing several periods. Conclusive evidence of the existence of a net directed current 



due to symmetric Levy flights is provide in Fig. 3. Consistent with 
current is observed. Note that, as in the case of Poissonian noise 



ig. 2(b), a net positive 
"up-hill" transport. 



i.e. transport the direction of the larger potential gradient, is observed. The dependence of 



the steady state current on the asymmetry parameter A for F = is illustrated in Fig. 4. 
Excellent agreement between the Langevin and the Fokker-Planck results is observed. For 
a fixed value of the Levy index a, the current is a monotonically decreasing function of the 
asymmetry and, as expected, for fixed A the current becomes smaller as a approaches the 
Gaussian limit a = 2. Consistent with the symmetry dxV{x,A) = —dxV{L — x, —A), the 
curves exhibit odd symmetry with respect the origin. Figure 5 shows the dependence of the 
steady state current on the Levy index a in the (3 = symmetric case. For large the 
current is a monotonically decreasing function of a. However, for small x the relation is 
non-monotonic and, depending on A, there is an optimal value of a that yields the maximum 
current. It is also instructive to explore the variation of the current with a fixing the scale 
factor cr = x^^" of the Levy distributions in Eq.Q. In this case the dependence seems 
to be always monotonic. A problem of interest in ratchet transport is the control of the 
magnitude and direction of the current. Figure 6 shows how this can be accomplished by 
manipulating the Levy noise asymmetry /3 and the potential asymmetry A. Consistent with 
the Langevin simulation in Fig. 2(c), in Fig. 6(a) there is a net current with A = and 
f3 = 0.25. Comparing with Fig. 3, for which A = —0.274 and /? = 0, it is observed that 
the noise asymmetric can mimic the potential asymmetry. As Fig. 6(b) shows, for a given 
ratchet potential, it is possible to find a compensating value of (3 for which the current 
vanishes. Moreover, as Fig. 6(c) shows, a further increase of the Levy noise skewness leads 
to a current reversal, a result consistent with the corresponding Langevin simulation shown 
in Fig. 2 (d). 

Summarizing. In this letter we have presented numerical evidence of directed transport 
driven by symmetric Levy noise in time-independent ratchet potentials in the absence of an 
external tilting force. In the limit a = 2, the noise becomes Gaussian, the fiuctuations are 
in thermodynamic equilibrium, and consistent with the second law of thermodynamics the 
current vanishes. However, for a ^ 2, the Levy noise drives the system out of thermodynamic 
equilibrium and an up-hill net current is generated. For small values of x and a fixed potential 
asymmetry, there is an optimal value of a yielding the maximum current. The direction 
and magnitude of the current can be manipulated by changing the Levy noise asymmetry 
and the potential asymmetry. As an application, we conjecture that a recently proposed 
ratchet pinch mechanism in magnetically confined plasmas might be present even in the case 
of static electrostatic fiuctuations. 
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FIGURE CAPTIONS 



FIG. 1. Effective ratchet potential Veff = V{x) — Fx for F = —1 according to Eq. (2), with 
A = -0.274 (dashed hne), A = 0.2 (dotted hne), and A = -0.5 (dashed-dotted 
hne). The sohd hne corresponds to the ratchet potential discussed in Ref. (see 
Fig. 2. 4) . The inset shows the well-known dependence on F in the a = 2 Gaussian 
case according to the Langevin simulations (dots) with A = —0.274, and the result 
reported in Ref. (solid line). 

FIG. 2. Typical trajectories according to the numerical solution of the overdamped Langevin 
Eq.([T]) with L = 1, X = 0.5 and: (a) a = 2 (Gaussian), A = -0.274; (b) a = 1.50, 
/3 = 0, A = -0.274; (c) a = 1.50, p = 0.25, A = 0; (d) a = 1.50 p = -0.5, 
A = -0.274. 

FIG. 3. Evidence of ratchet current in the Fractional Fokker-Planck model with a = 1.5, 
f3 = 0, X = 0.5, F = and A = —0.274. The solid line in the main figure shows 
the mean < x > with respect to x = which folows a linear scaling (dashed line) 
corresponding to a current < x >= 0.145. The inset at the top shows the final and 
initial (dashed line) probability distribution function superimposed with the ratchet 
potential in Eq. ([2]). The inset at the bottom shows the profile of the probability 
distribution function at the final time. 

FIG. 4. Steady current as a function of ratchet potential asymmetry A, for f3 = 0, x = 0.5 
and Levy indices a = 1.5, a = 1.75 and a = 1.90. The solid line with dots denote the 
results according to the Langevin model and the circles and crosses the Fractional 
Fokker-Planck results. The inset shows the dependence of the current on F for 
A = 0.2 and A = —0.2. The curve in the middle is the Gaussian case, which shows 
no dependence on the sign of A. The top (bottom) curve is the a = 1.5 Levy result 
for A = -0.2 {A = 0.2). 

FIG. 5. Steady state current < x > versus Levy index a for (a) x = 0.5, (b) x = 0.05, (c) 
a = 0.707, and (d) a = 0.0707. Curves 1, 2, 3, and 4 correspond to A = -0.274, 
A = —0.4, A = —0.5 and A = —0.6 respectively. 

FIG. 6. Dependence of ratchet current on Levy noise asymmetry in the fractional Fokker- 
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Planck equation. Panels (a), (b) and (c) show the probability distribution function 
for a = 1.5 and different values of /3, and A. The fourth panel shows the correspond- 
ing time evolution of the mean < x >. 
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